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We report here the frequency-dependent optical response of a nearly ferroelectric superconductor 
(NFE-SC), like Na x W03, or n-SrTi03. From uj = 0, up to a critical frequency uj c \, Meissner-like 
field extinction occurs as in a usual superconductor. In the "anomalous" region uj c i < uj < ujto 
(the TO frequency of the NFE) the field propagates as in a dielectric. For uj > ujto Meissner 
damping recurrs, then at higher uj, dielectriclike propagation is found. New optical anomalies 
(peaks in transmission) are predicted for a film of NFE-SC below the soft mode frequency ujto- 
These effects are due to including the contribution of the NFE displacement current in the solution 
of the coupled Maxwell, lattice and London equations. Quantitative calculations are given for the 
NFE-SC n-SrTi0 3 to illustrate this. 

PACS numbers: 74.20.Dc; 74.25. Nf; 78.20.Bh 



I. INTRODUCTION 

This paper predicts novel effects in the electromag- 
netic response of a material which exhibits supercon- 
ductivity (SC), and is in a " nearly- ferroelectric" (NFE) 
state. We use the description " nearly- ferroelectric" to 
mean a material which is a typical soft-mode phonon 
Cochran- Anderson system , with high-static dielectric 
constant. Materials in this class include the sodium tung- 
sten bronze Na^WOa and n- (or p-) doped SrTi03 sys- 
tems. 

Recently [2a] a sodium tungsten bronze Na K W03 with 
x ~ 0.05 was reported to be a high temperature super- 
conductor with T c ~ 90.fr. It is noteworthy that for 
0.1 < x < 1, Matthcis and collaborators long ago re- 
ported superconductivity but with T c ~ 3 — 5K [2b]. 
Also Mattheis and Wood [2b] first reported ferroelec- 
tricity in the host WO3 system. Another material, n- 
(or p-) doped SrTiOa exhibits superconductivity, with 
T c ~ 1 - 3K 3]. The host SrTi0 3 is known to be a 
" nearly- ferroelectric" material, where the static dielec- 
tric constant e(0) has been measured as ~ 10 4 in the 
temperature range where superconductivity occurs Q. 
In addition to these two materials, Weger and collabo- 
rators have recently proposed that in the high tempera- 
ture cuprate superconductors, a new multicritical point 
occurs due to an underlying nearly-ferroelectric instabil- 
ity which renormalizes the electron-phonon and electron- 
electron interaction and enhances T c 0, |(| . 

In the present work we will assume the "host" ionic lat- 
tice is a soft-mode NFE, with high dielectric coefficient 
e(0). To clarify the scope of our work we also assume 
that all the carriers injected in the host by doping are 
fully condensed into s-wave Cooper pairs in that lattice, 
so there are no "free" electrons. This assumption is con- 
sistent with earlier work on such systems by Cohen and 
Koonce :3J who presented a strong-coupling theory ap- 
plied to superconductivity in SrTi03 (see also Ref. |4{). 
Hence the free charge density p e = 0, and free current 



density J e = 0. 

In our present phenomenological approach, we solve 
the coupled equations for: a) the transverse electromag- 
netic field (Maxwell equations); b) soft-mode lattice vi- 
brations (Born-Huang equations); c) the superconduct- 
ing electrons in London (local) approximation. We re- 
port dispersive structure, when the electromagnetic fre- 
quency is tuned through the resonant frequency region 
near the basic ujto and ojlo frequencies of the soft mode. 
The theory is applied to a scmiinfinite medium, and to 
a film in order to calculate the reflectivity and trans- 
mittivity (impedance). We illustrate the results quan- 
titatively, using experimental data for superconducting 
doped SrTiOs, which is one of the materials for which all 
the needed lattice, electronic and superconducting data 
are available. 

We use London local electrodynamics for the super- 
conducting sector(SC), rather than a nonlocal BCS or 
Pippard approach. These NFE-SC materials are all type 
II (Xl > £, where Xt, is the London penetration depth, 
and £ is the coherence length) and in some cases they 
are strongly type II (Ax, >> £) [7]. It is well known that 
London electrodynamics is valid in these cases. We also 
do not consider in the present work effects due to vortex 
creation in these systems. 



II. MODEL AND DISPERSION 

At the outset we note that all the relevant dynamical 
fields considered, are transverse. These are the electro- 
magnetic fields, the lattice vibrations, and the London 
supercurrent. Maxwell's equations for the macroscopic 
fields are: 



V - B = 0, 



c ot 



V - D = 0, 



VxH-if = 0. (1) 

c at 
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with the constitutive equations of the medium taken as 

B = H + 4ttM and D = E + 4ttP = e(w)E. (2) 

For the host ionic lattice we assume a "diatomic" basis, 
with w the relative displacement vector of the (±) ions 
(w = u + — u_). The equations of motion of the lattice, 
omitting damping, are [8|: 



= 6nw + 6i 2 E, P = 6 2 iw + 6 22 E, (3) 



dt 2 



where bij are frequency-dependent coefficients and the 
electric field E and the polarization P are the same as 
in Eqns.(l) and (2). Now we seek a plane wave solu- 
tion proportional to exp(ik • r — iujt) for the Eqs. (3) in 
the long-wave approximation. We follow Huang, and for 
these transverse waves, in the ionic medium with a single 
resonance frequency ujto we obtain the usual expression 
for the dielectric function 



E ( W ) = 



{uj 2 



(4) 
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FIG. 1: Solutions of the dispersion equation (9) for zero 
damping. Frequency intervals uj c i < uj < ujto and uj > uj C 2 
correspond to real root of the Eq. (9). 
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where the expressions for the coefficients were given 
in Refs. 0. In media such as the perovskites, bronzes, 
or cuprates we have to take account of the other (non- 
soft) LO and TO modes by renormalizing to some "ef- 
fective" value e'^ 0, where e'^ = e ao W i uj 2 Li /ui Ti , (the 
primed product over "i" includes all oscillations except 
the soft mode) and we suppose the frequency uj is far 
below all frequencies ujm and wn except the soft mode. 
In what follows we use Eq. (4) with replaced by e^. 
The generalized Lyddane-Sachs- Teller (LST) relation is 
[10]: 



£q 



(5) 



Note that all oscillators, including i = (the soft mode) 
are included in Eq. (5) (no prime on the product). Fi- 
nally we take the London equation in the form: 

V x M = -J s = — tt-A = -—^—^A. (6) 

c c z A m*cr 

Here A is the vector potential for the electromagnetic 
field: 



VxA = B; -I^1=E. 

c dt 



(7) 



The London gauge has been used, and in Eq. (6) n s is the 
superconducting electron densityand m* is the electron 
effective mass. From Eqns. (1),(2),(4),(6) and (7) we 
obtain the equation for the electric field 



v- ,- ^ 1 , ^ 2 E 4tt 



c 2 A 



(8) 



The field equations describe the coupled electromagnetic 
field, the soft mode NFE oscillator, and the London su- 
percurrent. 

We now look for a transverse (T) plane-wave solution 
for E, and recalling that V x V x E = — V 2 E for such so- 
lutions, we obtain our new result, the dispersion equation 
for T waves 



-s{uj) 



(9) 



If the London term (last on the right-hand side) is ab- 
sent, then Eq. (9) is the usual equation for a transverse 
phonon polariton |8(; on the other hand, if the displace- 
ment term due to the polarization is absent (first term 
on the right-hand side) then the field attenuates in the 
distance A^ 1 = a/47t/c 2 A where Al is the London pen- 
etration depth Jjj. Owing to the dispersive form (4) for 
s(uj), interesting physics emerges as uj is increased from 
zero. Depending on the sign of the right-hand side of 
Eq. (9), k will either be real and electromagnetic waves 
will propagate, or imaginary so the wave will damp with 
frequency dependent penetration depth where 



Xl = X L /{\l-e{uj)\luj 2 /c 2 \) 



1/2 



(10) 



Figure 1 illustrates this. At uj = the usual London pen- 
etration depth \l is found. As the frequency increases 
\* l {uj) grows, up to a frequency uj c i, where k 2 (yj c \) = 
and \* l {uj c \) is infinite, so the field is uniform inside 
the medium. For uj c \ < uj < ujtq, k 2 > and the 
field propagates in the medium as in a dielectric. For 
ujto < ^ < ^JLQi k 2 < and again there is a fre- 
quency dependent A2(w), i.e. Meissner effect. Exactly 
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FIG. 2: Solutions of the dispersion equation (9) for zero 
damping. Frequency intervals ui c i < w < luto and ui > lu C 2 
correspond to real root of the Eq. (9). 



The roots of this equation are, 



^1,2 — <^TCM 

Here 
Ki = 1 



2r 2 



± -« N 2 



4r 2 H! 



\ 2 t4 (l + a 2 ) 2y 2 w 2 (l + a 2 ) 



(13) 



1 



1 + a 2 ' 



H 2 = 1 - 



1 



1 + a 2 ' 



% 2 = ^jrXcA 

47T 



is a dimensionless parameter which can take on values of 
the order of unity in nearly ferroelectric materials such 
as n-SrTi0 3 (see below). 

As we demonstrate below, when we calculate the ac- 
tual size of the region Slu — luto^^cI for n-SrTi03 using 
measured parameters including T, Slu is, we believe, suf- 
ficiently wide to enable measurement of our predicted 
effects in a film of NFE-SC. First we consider a semi- 
infinite medium, then the film. 



III. SEMIINFINITE MEDIUM 



at ujlo, £(<^lo) = 0, so X* l (ujlo) coincides with the Lon- 
don penetration depth. For lulq < uj < lu C 2, k 2 < 
giving increasing X^(lu), and at k 2 (uj C 2) = 0, a uniform 
field. Then, for u> > ui C 2, k 2 > and the wave propagates 
in the medium again, as in a dielectric. 

In summary, our dispersion equation (9) as a func- 
tion of frequency gives alternately regions of Meissner- 
like frequency-dependent penetration depth X* l (lu) that 
are superconducting with the magnetic field B excluded 
for x > \* L , then changing to regions of field propagation 
with w-dependent real wave number, as in a dielectric 
medium. This alternating dielectric and Meissner behav- 
ior suggests that our model has a zero-temperature phase 
transition between normal dielectric and superconductor 
phases driven by the electromagnetic frequency to, i.e. a 
type of " quantum phase transition" [12| . 

Up to this point damping was ignored . To examine 
this we write the dielectric function (4) allowing for a 
phcnomenological damping coefficient T (Ref. 



nng t< 

my- 



K - 



(oj^q — lu 2 ) — 2iujT 



(11) 



As shown on Fig. 2, including "moderate" damping 
r of the soft-mode oscillator changes the results quan- 
titatively but the above qualitative features remain the 
same. 

Our first step is to determine the extent of the impor- 
tant frequency interval lu c ± < lu < loto- The frequency 
LU c i is the smallest root of the Eq. (9) for k = 0, 



4tt 1 
= -t — J" 



(12) 



We consider a semiinfinite NFE-SC medium occupying 
the half-space z > 0. An incident electromagnetic wave 
propagates in the z direction with electric and magnetic 
fields E(r) and B(r) polarized along x and y axes respec- 
tively. The Maxwell boundary conditions at the interface 
plane z = are [14j 

E lx {0) + E rx (0) - E tx (0) = 0; 

4-7T 

E ix (0 )-E rx {0 )-— E tx (0 =0. (14) 
cZ 

Here £^(0), -E r:c (0), E tx (0) are the incident, reflected and 
transmitted electric fields, respectively and Z is the sur- 
face impedance. First consider the case without damping 
(r = 0), where we have 



z _ 4tt E tx (0) 
c B ty (0) 



dk 



±(\* L k) 2 



(15) 



We defined the effective penetration depth A£ by Eq. 
(10). In the expression (15) the upper signs are used for 
frequency ranges where the material exhibits Meissner- 
like behavior [k 2 < 0), while the lower signs apply to 
the frequency regime of dielectric- like behavior [k 2 > 0). 
After integrating over fc, we obtain 



Z= — 



A£ when lu < lo c \ or luto < lu < lu C 2\ 



and 



Z = — tt^l when lu c i < lu < luto, or lu > lu C 2- (16) 

Using these results we can solve the equations (14) and 
calculate the reflection coefficient for the Meissncr-likc 



4 



(k 2 < 0) or dielectric-like (fc 2 > 0). frequency ranges. 
We obtain the following: 



R = 



E rx {0) 



k X* LZ l 
+ 1 



k X) 



Ko 



(17) 



for the anomalous frequency ranges (u> c i < lu < luto\ u 
> ll>c2)', and R = 1 when lu < w c i or luto < u> < ^> c i- Of 
course the latter confirms a well-known result of electro- 
dynamics of superconductors, i.e., when we do not con- 
sider vortices, the electromagnetic field cannot penetrate 
into the volume of a superconducting medium. For the 
"anomalous" intervals, e.g., u> c i < lu < luto, the reflec- 
tion coefficient tends to unity near the limiting points lu c \ 
and loto because X* L tends to infinity when lu — ► lu c \ and 
A^ tends to zero when lu — ► luto- However in the interior 
of this frequency range R takes values which are consid- 
erably smaller than unity. This means that the electro- 
magnetic field can penetrate into a semiinfinite NFE-SC 
medium within the "anomalous" frequency range where 
the Meissner effect is suppressed. We did not include the 
effect of damping which will be discussed below. The 
damping effect is small in the practical case of n-SrTi(93. 



IV. THIN FILM OF NFE-SC 

We turn to the case of a slab or film of thickness L 
so that the material occupies the region < z < L. 
As previously, we take the electromagnetic field normally 
incident on the interface at z = 0, propagating in the z 
direction with the electric and magnetic fields depending 
only on z and parallel to x and y respectively. To proceed, 
we expand the electric field in the medium in a Fourier 
series, as done in Refs. 0, 



E(z) 



^ \ ~~ \ Sno ) En C0S ( kNZ ) 

N=Q ^ ' 

En sm(kiyz). 

N=l 



(18) 



Here 



En = I E(z) cos(fcjv z)dz; 
Iq 



TTiV 



En — / E(z) sin(kNz)dz; fc^r 

Jo L 

Using the dispersion equation (9), we get 



E 



N 



\*2 



gjL[E'(0)-(-l) N E>(L)] 



L n N 



En — 



l^^^[E(0)-(-1) n E(L)].(19) 

, N / 



The upper sign in the denominators of Eq. (19) is used 
for those frequency ranges where k 2 < 0, and the lower 
sign corresponds to the frequency ranges where k 2 > 0. 

Since the new results emerge in the region 5u> — loci — 
lotoi when k 2 > we consider this case in this Section 
[lower sign in (19)], and for the present we take T = 0. 
In this frequency range the material exhibits dielectriclike 
behavior. Substituting Eq. (19) into Eq. (18) we get 



£?x(0) = ^S,(0)oot(A 



iuj\ 



—By(L) sin" 1 



E X {L) = ^B y (0) sin" 1 ^y^B y (L) cot (Jl 

(20) 

The solution of this system of equations for the magnetic 
induction is: 

B „(o) = -^ E ,(o)cot(i)-^,( £)sta -'(A) i 



B , (L ) = _^ E ,(o )sta -.(A)__^ E , (i)cot ^ 

(21) 

The boundary conditions on the interfaces z = and 
z = L can be written as follows: 



E ix (0) + E rx (0) = E tx (0); 



E ix (0) - E rx (0) = B y (0); 



E X (L) = By(L). 



(22) 



Using these boundary conditions and expressions (21) for 
the magnitude of the magnetic field on the interfaces z = 
and z — L we arrive at the following expressions for 
the reflection and transmission coefficients: 



R = 



1 



T 



I + P 2 (lj) sin- 2 (L/\* L ) 

p 2 (u;)sm- 2 (L/\i) 
l + p 2 {L0)s\n- 2 {L/\iy 



where the the function p(u) has the form: 

, , 2ck* 1 
p(uj) = 

LU 1 — C Z K* Z /UJ Z 



(23) 
(24) 

(25) 



Here k* = l/X* L is the solution of the dispersion equa- 
tion (9) in the "anomalous" frequency range loci < 
lu < loto- For frequencies very close to luto we have 
k = {u>/c)yfe which corresponds to a polariton-phonon 
transverse wave. 

We now note an interesting result following from Eqs. 
(23) and (24), when the effective London penetration 
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FIG. 3: Solutions of the dispersion equation (9) for zero 
damping. Frequency intervals ui c i < lo < loto and w > uj C 2 
correspond to real root of the Eq. (9) . 



depth X* L , and the Fourier component k = kjy are in reso- 
nance. That is, for X* L = L/(nN) or k = k^, sin(L/A*) in 
Eqs. (23) and (24) equals zero and we get for the reflec- 
tion and transmission coefficients the values R — 0; T — 
1, which means total transparency of the slab at A: at (See 
Fig. 3). 

Total transparency can also occur at a frequency for 
which k = k = u/c. At this frequency we have 
e(oj) — 4it/oj 2 A = 1 so it is located close to the left 
boundary frequency u) c \. When ko = k the function p(u) 
tends to infinity and again we get R = and T = 1. 



To explain this we can treat the quantity yj e{uS) = 
\J e{uo) — A-n/uj 2 K = n(uj) as an effective reflective index 
of the slab material. When it equals unity the electro- 
magnetic field from outside passes through the slab with- 
out reflection. 

Our results (23) and (24) were obtained, neglecting 
damping. For a nonzero damping constant T we have to 
replace our expressions (23) and (24) near the resonances 
by the following: 



R 



1 + p 2 (uj) + S 2 (u) sin- 2 {L/Xl'y 

p 2 (u J )sm- 2 (L/\j') 
1 + p 2 (co) + S 2 (lo) sin- 2 (i/A* ') ' 



(26) 
(27) 



S(u) - p(w) cos (J^J + v/T+TMsinh (J^jj 

(28) 

andl/Ai = l/Ai' + i/Ai". 

We can use the following approximation for X* L ' and 
XI" : 




LU 2 1 



2loT 



UJ^ — UJ< 



1/2 



1 



(29) 
(30) 



where 



It was shown before that damping in the lattice sys- 
tem reduces the anomalous frequency range Su>. The fre- 
quency interval w c i < uo < ujto has to be replaced by the 
narrower interval cj\ < u < uu 2 - The imaginary part of 
1/A^ depends on frequency and increases when the latter 
increases. As we will illustrate below, moderate damp- 
ing, as in n-SrTiC>3 will only slightly modify the sharp 
resonances which occur at k = k^ and k — fcrj. 

To the best of our knowledge the general behavior we 
predict for an NFE-SC system has not been observed. In 
order to examine the practicality of making such obser- 
vations we have calculated the effects quantitatively for 
n-doped SrTiC>3. This material was much studied, as a 
nearly-ferroelectric, and as a superconductor, prior to the 
era of high T c , so that the needed physical parameters 
are available. 



V. APPLICATION TO N-SRTIO3, A 
PROTOTYPE NFE-SC 

In order to evaluate the magnitudes and locations of 
the predicted effects, we require lattice, dielectric, and 
electronic data. So far as we can determine such mea- 
sured data is only available for n-SrTiC>3, and we will use 
it below. Nearly ferroelectric SrTi0 3 has been studied for 
some time. The temperature dependent static dielectric 
coefficient e(0) has been measured in a wide temperature 
range, extending down to about IK [4]. The frequency 
dependent is known at various temperatures [4]. Lat- 
tice dynamic studies were carried out based on the shell 
model [13], and neutron scattering [13], and infra-red and 
Raman (including electric field dependence) optical prop- 
erties have been measured [9] . As a result the lattice dy- 
namical parameters are known, and the generalized LST 
relation (5) has been verified. From the line-widths in 
infra-red and Raman properties [9] , the damping param- 
eter T for the soft mode lattice frequency loto is known. 
Hence all the physical parameters of the dielectric host 
are actually known. The temperature independent quan- 
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tities are 

Wui 2 

£00=5.5, uj LO = 5.2 x 10 12 s"\ -4^=4.1. 

Using the LST relation we have at ~1K 
ujto = 1-6 x lO 1 ^- 1 

and at IK 

e(0) - 2.25 x 10 4 . 

The damping constant for the soft TO mode as in eqn(13) 
is measured from optical line width [9|: 

r 2 /u>T ~ °-2 

at the temperature of interest. (Actually the measure- 
ment was at a higher temperature than the region of T c 
so the relevant T for our work will be smaller, but we use 
the above value.) 

For the needed electronic parameters we use the effec- 
tive mass m* ~ 10m e where m e is the free electron mass. 
A range of n doping from 10 17 to 10 21 cm -3 was used 
for SC n-SrTiC>3. We choose n s =9x 10 17 cm~ 3 for our 
calculation. We recall that n-doped SrTiC>3 exhibits con- 
centration dependent T c , with T c = 0.3K in the optimum 
doped system |3j. In this temperature range e(0) has the 
high value given above, as was noted also by Cohen 

We now determine the anomalous frequency interval 
just below ujto- We first need the parameter a 2 = 
uIqe'^K/A-k [see Eq. (13)]. For SrTi0 3 with n s = 
9 x 10 17 cm~ 3 , and to* w 10m e , we get a 2 w 2.1. Tak- 
ing, first, r = 0, and solving Eq. (12) we find ui c \ = 
0.9 x 10 n c _1 . This gives 5u) = uj to -uj c i = 0.43u; TO . Us- 
ing the measured T for SrTiC>3, we obtain that 0J2 shifts 
negligibly while lo c \ shifts slightly to ui' cl — 0.96x 10 11 c _1 , 
and the anomalous frequency region becomes slightly 
smaller 

Suj' = ujto — w cl ~ OAujto- 

We expect that this interval Suj' is still sufficiently wide 
to examine optical response experimentally by reflec- 
tion/transmission studies in n-SrTiC>3. The small effect 
of damping makes it likely also that T = is a good 
approximation for the calculations that follow. 

We now turn to calculating the reflection/transmission 
coefficients for a n-SrTiC>3 film of NFE-SC of thickness L, 
first for T = 0, using the formulas (23)-(25). In this case 
the results are shown on Fig. 3. We predict a "comb-like" 
structure of narrow perfect transmission T p = 1 spikes al- 
ternating with regions where T p — 0. The spikes occur at 
the resonance condition where A^ = L/ttN. Next we es- 
timate the effect including damping. For weak damping 
T 2 /uj to = 0.2 the value of the factor [2ujT /(uj tq ~uj 2 )] 1/2 
changes from 0.9 to 1.5 when u> runs over the interval be- 
tween u>i and ll>2- Thus, when the thickness of our film 



L is small or of the same order as the lu = London 
penetration depth Al (here Al ~ 1.8 x 10 -3 cm) mod- 
erately weak damping as in SrTi03 cannot significantly 
change the height of the peak in the transmission coef- 
ficient. Peaks arise due to the propagation of the trans- 
verse mode corresponding to the solution of the disper- 
sion equation (9). We believe that these peaks can be 
observed in any NFE-SC material in the appropriate fre- 
quency ranges, and at temperatures below T c . Such an 
observation would give support for theory proposed here. 

It is not possible to make realistic estimates of the 
interval Su) (or Suj' with damping) for Na^WOs, or for 
cuprates pending the availability of the necessary mate- 
rial constants as we had in n-SrTi03. The theory pre- 
sented will apply to any such NFE-SC material, so the 
experimental search for the resonance comb-like structure 
in T(w) just below ui = ujto will be highly valuable. 



VI. SUMMARY 

We solved the coupled "soft-mode" lattice dynamics 
(long wavelength) equations of motion, together with 
Maxwell equations for the associated electromagnetic 
field and the London equations for the supercurrent, in 
a NFE-SC. The analysis can apply to several material 
systems of current interest: sodium tungsten bronzes 
Na^WOs, doped n-SrTi03, and possibly high temper- 
ature cuprates. We note that several authors have dis- 
cussed microscopic theories for the superconductivity in 
n-SrTi03 (Ref. [15J) in the strong-coupling framework. 
However, our macroscopic approach has not previously 
been reported, to our knowledge. We assume all the free 
carriers are condensed in Cooper pairs, as in some earlier 
work on the SrTi0 3 (See Regs. 0, 0|). The resulting 
coupled modes can be considered as phonon-polaritons 
dressed by the supercurrent or electromagnetic waves in 
a London superconductor dressed by the transverse optic 
waves (TO phonons). 

A result, which we illustrated, is the alternation of 
the system response between Meissner-like superconduc- 
tor and a dielectric-like medium, as the incident electro- 
magnetic frequency is continuously varied. Particularly 
of interest is the "comb- like" series of transmission res- 
onance peaks we predict, which should be measurable. 
Such measurements would test our theory. One NFE-SC 
material of choice would be n-SrTi03 for which we made 
quantitative estimations. 
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APPENDIX: PARTITION OF ENERGY 

BETWEEN RADIATION, LATTICE 
POLARIZATION, AND LONDON FIELD 



Since all waves are transverse take E = (£^,0,0); B = 
(0, B v , 0); k = (0, 0, fe) giving B y = (ck/uj)E x . We elimi- 
nate H from (A.l), and combining all terms, we get 



It is useful to obtain the frequency-dependent parti- 
tion of energy between the three propagating fields. This 
can be seen as extending Huang's early treatment of the 
radiation-lattice coupling (phonon-polariton y|) to in- 
clude the London supercurrent, or conversely extending 
London treatment [16( of the energy momentum theorem, 
to include the dielectric polarization. 

Within our treatment all fields are transverse and we 
consider the infinite medium. We will obtain the conti- 
nuity equation for energy density U and Poynting energy 
flux vector s. From Maxwell's equations (1) we have: 



V-s 



1 dB 1 dB , „ , 

— H — -— E — , (A.l) 

4tt dt 4tt dt v ' 



where 



4?r 



[E x H]. 



(A.2) 



In order to simplify the right hand side of (Al) we use 
the lattice equation of motion (3) to obtain: 



f 



dB 



f 







2tt d ( d 



2c °° dt 



c dt\dt 



(A3) 

Then, using the first of equations (2), and equations (6) 
and (7) we have 



[V x Ml 



f 



1 



(A4) 



and 



dU 



1 



dB 1 dD 



1 
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dt 4tt dt 4tt dt 8tt 
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'dt 



E 



dB 1 

dt 



(A.10) 



It is simple to identify the various terms in the energy 
density U as belonging to radiation field, lattice polar- 
ization, and London supercurrent. Thus 



Urad = ^-£ooE 2 + -^B 2 , 



U 



LATTICE 



1,9 1 9 2 

-w + -w TO w , 



Use = 



1 1 

S^Ap 2 



(All) 



(A.12) 



(A.13) 



where fco = w/c. All these contributions to the energy 
density can be expressed in terms of E 2 in our plane- wave 
case, giving 



U LATTICE 



(gp - £oo)^r (^ 2 + Mjo) F 2 

87r(w 2 -w 2 ) 2 



(A.14) 



(A.15) 



[V x [V x Ml] = -V 2 M = — — B 

c 2 A 



(A. 5) an( i, as before: 



Now we will assume plane monochromatic electromag- 
netic wave propagation (~ exp(ik • r — iuit)). So 

d 

V 2 ~k 2 - — -iu. 

dt 



Then 



fc2M = -^ B < 



(A.6) 



Use = 



1 1 



(A.16) 



Before proceeding we note that the expression for the 
lattice term in energy density is unchanged from that 
found by Huang || for the phonon-polariton. Formally 
that is also true for the radiative term, however hidden 
in Urad is the effect of the London term via the changed 
dispersion equation (9), which includes A. Note that sub- 
stituting the dispersion equation (9) into Urad we get 



B = H 



or with B = fiH. we get: 
1 

- = 1 
A* 



and 



B. 



A|fc 2 



1 dB_ iig^B j_ 

c dt c n dt 2c 



(A.7) 



(A.8) 



A 2 k 2 J dt 



Ib 2 . 



(A.9) 



Urad = ^ ( £oc + e(w) - ™ ) E 2 (A.17) 
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and 
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(A.18) 
(A.19) 
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We may see that the energy-density contribution Use 
will exactly cancel the London contribution to the radi- 
ation energy density Urad- Thus, if we calculate the 
lattice (polarization) fraction of the energy we obtain ex- 
actly the result of Huang [j| despite the coupling of the 
London supercurrent. We have 



Plattice 



U LATTICE 



U LATTICE 
V TOTAL V LATTICE + Urad + Use 



1 



(uj 2 



TO 



-uj 2 )(uj 2 



LO 



J TO 



2 (uj 



2 

TO 



,2\2 



- U), 



TO\ w LO 



- ur 



(A20) 



This result shows that the mode described by the dis- 
persion equation (9) is a mixed photon-phonon mode 
within all frequency range w c \ < u> < ujto- We empha- 
size this result is only meaningful in the intervals where 
k 2 > 0. The lattice contribution to the energy enhances 
when the frequency increases and gets its maximum when 
ui goes to loto- At the frequencies close to ujto p ~ 1, so 
all the energy density is concentrated in the ion system. 
Taking into account that the frequency ujto is a bound- 
ary frequency for the "anomalous" frequency range, we 
can expect that the ion oscillations at the frequency ujto 
being excited by the electromagnetic field which can pen- 
etrate to a NFE-SC material within the " anomalous" fre- 
quency interval, will exist at frequencies larger than ujto 
where the material behaves as a superconductor. 
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